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Abstract. An optimization program is used to re-adjust the initial conditions, in order to reproduce as closely 
as possible the predictions of a complete ephemeris by using simplified equations of motion in the numerical 
integration. The adjustment of the initial conditions is illustrated in the transition from the DE406 complete 
long-range ephemeris to a Newtonian model considering only the Sun and the four major planets. It is also used 
to best reproduce this same DE406 ephemeris, based on post-Newtonian equations for a system of mass points 
and including the Moon and asteroids, by using a Newtonian calculation corrected by the Schwarzschild effects 
of the Sun and restricted to the ten major bodies of the solar system. 
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1. Introduction 

When integrating the equations of motion in classical or 
relativistic celestial mechanics, one has to know the val¬ 
ues of the parameters, in particular the masses and the 
initial conditions. The masses of the celestial bodies are 
obtained by combining different methods, including the 
analysis of the perturbations that they cause to the mo¬ 
tion of other bodies, and the fitting of the observed tra¬ 
jectory of a spacecraft in close approach. In the construc¬ 
tion of modern ephemerides, the initial conditions are ob¬ 
tained by le ast-squares fittings to large sets o f observa¬ 
tional data ( Newhall et al. 1983 ; Ficnga 1999| ). It is im¬ 
portant to realize that the optimal parameters, that lead 
to the minimum (least-squares) residual with respect to a 
set of observational data, depend on the precise model that 
is used. Thus, in the case that an analytical ephemeris is 
adjusted to a numerically-integrated one, the initial con¬ 
ditions taken from the numerical ephemeris have to be 


slightly modified (Arminjon 2000). A program for the ad¬ 
justment of the masses and the initial conditions has been 
built by the author in the latter context. It has been 
tested by adjusting a purely Newtonian calculation to the 
numerically-integrated ephemeris DE403 ([Standish et al. 


iterative minimization of the mean quadratic error, which 
needs to calculate the partial derivatives of the theoretical 
predictions with respect to the parameters. Several stan¬ 
dard methods may be used to compute these derivatives. 
In this program, they are calculated by finite differences: 
this involves making loops on the numerical solution of 
the equations of motion. 

In the present paper, this optimization program is used 
to re-adjust the initial conditions in the transition from 
the complete ’’long” ephemeris DE406 (Standish 1998) 


modified (Lestrade & Bretagnon 1982; Moisson 2000). If 
one aims at testing an alternative theory of gravity, one 
may expect that even the masses will have to be very 


1995), that is based on general-relativistic equations of 


to simplified models that are more tractable for integra¬ 
tions over very long times. Thus, the input data for the 
adjustment are taken from the complete ephemeris, and 
two simplified models of the solar system will be consid¬ 
ered: in Section 2, we present the adjustment of a model 
based on purely Newtonian equations, and limited to the 
Sun and the four major planets. In Section 3, the ini¬ 
tial conditions are re-adjusted for a model that consid¬ 
ers the nine major planets and accounts, moreover, for 
those general-relativistic effects that are caused by the Sun 
alone. For these two adjustments, the sets of the optimal 
initial conditions are provided for potential users, and the 
differences with the predictions of the reference ephemeris 
are illustrated. Moreover, the effect of extrapolating the 
calculation outside the fitting interval can be seen, be¬ 
cause the fitting intervals for the different planets are all 
significantly smaller than the 60 centuries of the DE406 
ephemeris. 


motion. The program is based on the Gauss algorithm for 
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Fig. 1. Comparison between DE406 and an adjusted Newtonian calculation restricted to the Sun and the four giant 
planets, over 60 centuries. Plotted are the differences between the heliocentric position vectors, distances, latitudes, 
and longitudes, in the sense Newtonian calculation minus DE406 ephenreris. Fig. la: Jupiter. 


2. Initial conditions for Newtonian integrations 
with the four major planets 

The initial conditions, at Julian Day 2451600.5 (26 
February 2000, 0H00), of a Newtonian calculation con¬ 
sidering five mass points: the Sun and the four giant plan¬ 
ets, were adjusted. The target time interval for the ad¬ 
justment was from -1000 to +2000. The input data (41 
heliocentric position vectors for each planet, equally dis¬ 
tributed in time for a given planet) were taken from the 
DE406 ephemeris, and spanned a good part of this inter¬ 
val: 1073 years for Jupiter, 2672 for Saturn, and 2913 years 
for Uranus and Neptune. The initial conditions were ini¬ 
tialized at the DE406 values, then adjusted, whereas the 
GM products were fixed at the DE406 values. The stan¬ 
dard deviation, between the input data of DE406 and the 
corresponding values obtained with the Newtonian calcu¬ 
lation restricted to the Sun and the four major planets, 
was reduced by a factor of 434 by the adjustment program. 
This ratio demonstrates the necessity of readjusting the 


initial conditions, as was done. Table [I] shows the values 
found for the initial conditions after this adjustment. 

Figure [[] shows the comparison between DE406 and 
this adjusted Newtonian calculation, over the 60 centuries 
(-3000 to +3000) covered by DE406. Heliocentric coordi¬ 
nates are being used for the figures and for the tables as 
well. The reference plane for the tables is the one which 
was used in the integrations, namely the Earth’s J2000 
equator. For the figures, the coordinate system has been 
rotated so that the reference plane is the invariable plane 
of the solar system (normal to the total angular momen¬ 
tum of the system). It can be seen that, using the adjusted 
initial conditions, the difference with DE406 remains very 
small, even over the extrapolated time range of 60 cen¬ 
turies. This is especially true if one remembers that the 
model used here is very strongly simplified as compared 
with the complete calculation of DE406. However, one ob¬ 
serves a sensible deterioration outside the target interval 
of the adjustment (which, we recall, was from —1000 to 
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Fig.1. b: Saturn. 


+2000), especially for Neptune. The most favourable be¬ 
havior in the extrapolation is found for Jupiter, for which 
the errors in distance and latitude increase merely linearly 
with time, whereas the longitude shows hardly any long- 
range drift. It is, of course, not possible to predict how 
the differences would evolve outside the 60 centuries of 
DE406. 


3. Initial conditions for calculations with the nine 
major planets, accounting for relativistic effects 
in the Schwarzschild field of the Sun 

The ephemerides series DExxx, constructed at the Jet 
Propulsion Laboratory, are based on Einstein-Infeld- 
Hoffmann-type equations of motion, and thus include 
general-relativistic effects at the first post-Newtonian ap¬ 
proximation in the solar system, the latter being con¬ 


sidered as a general system of point masses (Newhall el 


to the Newtonian acceleration of any given planet, that 
correction which is obtained when considering this planet 
as a test particle in the Schwarzschild field of the Sun [cf. 
Brumberg (1991)], thus writing the heliocentric accelera¬ 
tion of a planet as: 


du 

At 


= heliocentric acceleration in the Newtonian 


TV-body problem 


GM 


+ — (u.n)u ) 

c z 


GM 


3 / \ 2 : 

-(u.n) u 


(1) 


where G is the gravitation constant and c is the veloc¬ 
ity of light, M is the mass of the Sun, r = ||x|| = (x.x) 1 / 2 , 
n = x/||x||, and u = dx/df, x being the heliocentric po¬ 
sition vector. (Equation (l!) applies when Schwarzschild’s 


al. 1983). If one wishes to account merely for the major 


general-relativistic effects, one may consider instead the 
Newtonian TV-body problem for point masses, and add, 


solution is expressed in standard coordinates.) Lestrade 
fe Bretagnon (1982|) proposed an analytical integration of 


Eq. (|3j) , based on a Taylor expansion of Gauss’ perturba¬ 
tion equations for the osculating elements. They presented 
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Fig. 1. c: Uranus. 


numerical results for Mercury, and found very small differ¬ 
ences with the DEI02 ephemeris, after the integration con¬ 
stants of their integration had been adjusted on DEI02. 
In the present work, a numerical integration of Eq. ( 0 ) has 
been used, based on the Matlab solver ODE113. [For all 
calculations in this paper, the same numerical tolerances 
have been used in the ODE113 routine as those found opti¬ 
mal in our first test of this program, i.e. RclTol = 5 x 10“ 13 
and AbsTol = 10 -15 . Moreover, using the osculating ele¬ 
ments, some analytical perturbations of the Moon on the 
motion of the EMB (Bretagnon 198C) were taken into ac¬ 
count in the program.] We found only a marginal increase 
in the computation time when passing from the Newtonian 
heliocentric equations of motion to Eq. (jlj). The reasons 
are that: i) the program spends the greatest part of the 
time in calculating the Newtonian accelerations due to 
the other planets, and ii) when the ” Schwarzschild correc¬ 
tions” are added, the set of the step sizes automatically 
adopted by the solver in order to reach the given toler¬ 
ances, does not have a smaller average than has the set of 
the step sizes adopted with purely Newtonian equations 
(and with the same tolerances). The latter point means 


that the Schwarzschild corrections do not ’’stiffen” the dif¬ 
ferential system. 

However, the integration time is roughly 200 times 
greater with ten bodies including Mercury, than for five 
bodies with Jupiter as the quickest planet (the latter case 
was investigated in Section 2), in addition there are more 
initial conditions to adjust. For this reason, only ’’short- 
range” adjustments may be run easily with common com¬ 
puters. For the present adjustment, 21 positions per planet 
were taken from the DE406 ephemeris. As for the previ¬ 
ous adjustment, the time intervals spanned by these input 
data increased with the period of the planet: 4000 days 
for Mercury, and 64000 days (175 years) for the five outer¬ 
most planets. The adjustment program reduced the stan¬ 
dard deviation between the input data of DE406 and the 
corresponding values obtained with the calculation based 
on Eq. i), by a factor of 492. This again demonstrates 
the necessity of readjusting the initial conditions. Table || 
gives the initial conditions found after this adjustment. 

Figure shows the comparison between DE406 and 
this adjusted Newtonian calculation with Schwarzschild 
corrections, over the 60 centuries of DE406. If one first 









































































































































































M. Arminjon: Proper initial conditions for long-term integrations of the solar system 


5 


X 10“ 4 



Time (years) 



Fig. 1. d: Neptune. 

compares Figures 2e-f-g-h with Figures la-b-c-d, one ob¬ 
serves that, for the giant planets, the present calculation 
gives gains of accuracy (the latter being measured, e.g., by 
the norm SR of the radius vector difference) by a factor 
of: 20, 250, 300, and 3000, for Jupiter, Saturn, Uranus, 
and Neptune, respectively, as compared with the simpler 
calculation of Section 2. The poorest accuracy is found 
for Mercury (Fig. 2a). Over one century, the error is very 
small for Mercury also: 6 x 10 -8 au for SR and 0.03" for 
the longitude, but it increases more quickly with time for 
Mercury. This may be due, at least for one part, to the 
numerical integration error. 

The effect of changing the weights was considered: 
by taking weights that are proportional to the inverse 
of the data uncertainty, one reduces the standard devi¬ 
ations between input data and their recalculated values, 
by a factor two for Venus and the EMB, by 50% for Mars, 
and by 25% for Mercury, although the residuals for outer 
planets are all deteriorated (by factors of ten for Jupiter, 
and of several hundreds for Uranus, Neptune and Pluto). 
In order to still improve the agreement with a complete 
ephemeris, one could also consider other fine points, e.g. 


x 10' 5 



Time (years) 


Neptune, adjusted Newton vs DE406, 60 cy 



Time (years) 


whether the same form of metric (isotropic or harmonic, 
etc.) is being used in the Schwarzschild solution as in the 
post-Newtonian equations of motion for a system of mass 
points. But, if the accuracy over long time ranges is aimed 
at, one should first check the effect of improving the ac¬ 
curacy of the numerical integration (e.g. by passing from 
double to quadruple precision). 

4. Conclusion 

The adequate values of all parameters in the equations of 
motion depend on the precise form adopted for the latter 
equations. In particular, it is very important to re-adjust 
the initial conditions when one passes from a complete 
form of the equations of motion, as it is used to build 
reference ephemerides, to a simplified form, which may 
be more adequate for long-range calculations. Although it 
may be too time-consuming to optimize the initial con¬ 
ditions so as to get the smallest residual for the whole 
time-range to be considered, one can get ’’reasonably ex- 
trapolable” initial conditions, by optimizing them for a 
time interval that contains at least a few periods of even 
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Fig. 2. Comparison between DE406 and an adjusted calculation restricted to the Sun and the nine major planets, 
and including Schwarzschild corrections of the Sun, over 60 centuries. Plotted are the differences between the helio¬ 
centric position vectors, distances, latitudes, and longitudes, in the sense (Newtonian calculation with Schwarzschild 
corrections) minus (DE406 ephemeris). Fig. 2a: Mercury. 


the slowest body considered. Provided one thus re-adjusts 
the initial conditions, it is possible to reproduce quite ac¬ 
curately the corrections of post-Newtonian general relativ¬ 
ity, by adding the Schwarzschild corrections, that account 
for the post-Newtonian effects of the Sun alone. This does 
not significantly increase the computer time as compared 
with the purely Newtonian equations of motion. A very 
precise comparison between the complete and simplified 
ways of describing post-Newtonian effects in the solar sys¬ 
tem could be obtained only if the two models were iden¬ 
tical in the other respects, which was not the case here. 
Including the masses in the solved-for parameters brought 
only marginal improvements in the models for which it was 
considered. 
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Fig. 2. c: Earth-Moon Barycenter. 


Table 1. Initial conditions at JD 2451600.5 found for a Newtonian model containing only the Sun and the 4 giant 
planets after the adjustment (heliocentric equatorial Cartesian coordinates). 


Positions (au) 
Jupiter 
Saturn 
Uranus 
Neptune 

+3.733061134680E+00 
+6.164409988928E+00 
+1.457963610839E+01 
+1.695475862942E+01 

+3.052435941481E+00 

+6.366764207406E+00 

-1.236887103144E+01 

-2.288718303371E+01 

+1.217433227240E+00 

+2.364527399872E+00 

-5.623600018369E+00 

-9.789935831446E+00 

Velocities (au/d) 
Jupiter 
Saturn 
Uranus 
Neptune 

-5.086558851107E-03 

-4.426842485459E-03 

+2.647476590726E-03 

+2.568634013660E-03 

+5.493661581130E-03 

+3.394095269623E-03 

+2.487495861475E-03 

+1.681857805109E-03 

+2.478691066898E-03 

+1.592276937224E-03 

+1.052019949681E-03 

+6.245754662095E-04 
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Table 2. Initial conditions at JD 2 451600.5 found for the model based on 


Eq. (0). 


Positions (au) 
Mercury 
Venus 

Earth-Moon Baryc. 
Mars 
Jupiter 
Saturn 
Uranus 
Neptune 
Pluto 


-2.503321047836E-01 
+1.747780055994E-02 
-9.091916173950E-01 
+ 1.203018828754E+00 
+3.733076999471E+00 
+6.164433062913E+00 
+1.457964661868E+01 
+1.695491139909E+01 
-9.707098450131E+00 


+1.873217481656E-01 

-6.624210296743E-01 

+3.592925969244E-01 

+7.270712989688E-01 

+3.052424824299E+00 

+6.366775402981E+00 

-1.236891078519E+01 

-2.288713988623E+01 

-2.804098175319E+01 


+1.260230112145E-01 

-2.991203277122E-01 

+1.557729610506E-01 

+3.009561427569E-01 

+1.217426663570E+00 

+2.364531109847E+00 

-5.623617280033E+00 

-9.789921035251E+00 

-5.823808919246E+00 


Velocities (au/d) 
Mercury 
Venus 

Earth-Moon Baryc. 
Mars 
Jupiter 
Saturn 
Uranus 
Neptune 
Pluto 


-2.438808424736E-02 
+2.008547034175E-02 
-7.085843239142E-03 
-7.124453943885E-03 
-5.086540617947E-03 
-4.426823593779E-03 
+2.647505630327E-03 
+2.568651772461E-03 
+3.034112963576E-03 


-1.850224608274E-02 
+8.365454832702E-04 
-1.455634327653E-02 
+1.166307407692E-02 
+5.493643783389E-03 
+3.394060157503E-03 
+2.487457379099E-03 
+1.681832388267E-03 
-1.111317562971E-03 


-7.353811537540E-03 
-8.947888514893E-04 
-6.310912842359E-03 
+5.542098698449E-03 
+2.478685100749E-03 
+1.592261423092E-03 
+1.052000252243E-03 
+6.245613982833E-04 
-1.261841468083E-03 
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Mars, "Newton+Schwarzschild" vs DE406, 60 cy 



Fig. 2. d: Mars. 
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